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SUMMARY 


The  finite  difference  based  method  of  characteristics  model  for  unsteady 
partially  filled  pipe  flow  was  extended  to  include  the  model  for  the  stack 
to  horizontal  drain  entry  boundary  condition.  The  conditions  at  drain  entry 
are  defined  in  terms  of  the  energy  of  the  terminal  annular  flow  velocity  in 
the  stack,  together  with  an  appropriate  loss  coefficient  function  at  the  entry. 
The  hydraulic  solutions  link  the  branch  drains,  fittings,  vertical  soil  stack, 
and  building  drain.  The  analysis  permits  any  combination  of  drainage  load 
patterns  from  simultaneous,  overlapping,  or  sequence  of  discharge  events. 
Preliminary  simulations  utilizing  this  model  indicate  that  the  modeling  tech- 
nique extends  the  existing  horizontal  network  analysis  program  for  determination 
of  multistory  building  drainage  systems.  The  sizing  procedure  provides  the 
hydraulic  capacity  of  drains  for  specified  pipe  sizes,  pipe  pitch,  and  wall 
roughness  factors. 

Keywords:  Drain  stack  to  horizontal  drain  flows;  drain  system  design; 

plumbing  drainage;  unsteady  flow  in  drains. 
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SUMMARY 


The  finite  difference  based  method  of  characteristics  model  for  unsteady 
partially  filled  pipe  flow  was  extended  to  include  the  model  for  the  stack 
to  horizontal  drain  entry  boundary  condition.  The  conditions  at  drain  entry 
are  defined  in  terms  of  the  energy  of  the  terminal  annular  flow  velocity  in 
the  stack,  together  with  an  appropriate  loss  coefficient  function  at  the  entry. 
The  hydraulic  solutions  link  the  branch  drains,  fittings,  vertical  soil  stack, 
and  building  drain.  The  analysis  permits  any  combination  of  drainage  load 
patterns  from  simultaneous,  overlapping,  or  sequence  of  discharge  events. 
Preliminary  simulations  utilizing  this  model  indicate  that  the  modeling  tech- 
nique extends  the  existing  horizontal  network  analysis  program  for  determination 
of  multistory  building  drainage  systems.  The  sizing  procedure  provides  the 
hydraulic  capacity  of  drains  for  specified  pipe  sizes,  pipe  pitch,  and  wall 
roughness  factors. 
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1 . INTRODUCTION 


The  unsteady  flow  equations  defining  partially  filled  pipe  flow  may  be  solved 
by  means  of  the  finite  difference  method  of  characteristics.  Programs  to 
utilize  this  technique  have  been  developed  and  fully  described  in  previous 
progress  reports  and  published  papers  [1]. 


The  development  of  the  basic  equations  and  their  application  in  finite 
difference  form  to  describe  unsteady  partially  filled  pipe  flow  are  illustrated 
by  figures  1 and  2.  In  common  with  other  finite  difference  applications  the 
methods  require  a knowledge  of  both  initial  conditions  and  the  system  boundary 
equations  to  allow  the  solution  to  proceed.  For  a single  pipe  as  shown  in 
figure  2,  under  supercritical  flow  conditions,  the  exit  boundary  is  supplied 
by  the  C+C~  characteristics  from  the  preceding  upstream  section.  In  addition 
to  the  boundary  conditions  the  finite  difference  solution  also  requires  know- 
ledge of  the  initial  flow  depth  and  velocity  along  the  pipe  length  at  the 
start  of  the  simulation.  Normally  a small  steady  flow  is  assumed  along  the 
whole  pipe  length  prior  to  the  introduction  of  the  drain  load  profile. 


The  flow  profiles  entering  the  drain  may  be  assumed  for  initial  case  studies. 
For  example,  a high  energy  steady  inflow  will  become  progressively  deeper 
along  the  pipe's  length  as  its  energy  is  dissipated.  The  flow  depth  will 
eventually  reach  the  normal  steady  uniform  flow  depth,  however  the  distance 
required  for  this  will  depend  on  the  entry  energy,  the  flow  rate  and  pipe 
parameters  such  as  diameter,  roughness,  slope,  etc.  The  entry  condition  at 
the  head  of  the  drain  depends  on  the  inflow  conditions  resulting  from  the  flow 
upstream  of  the  entry  section.  In  this  report  the  upstream  flow  conditions 
are  derived  to  represent  fully  developed  vertical  stack  flow  entering  the 
"horizontal"  drain.  slope  etc.  Determination  of  this  initial  flow  depth 
profile  as  a function  of  time  at  the  horizontal  drain  entry  will  be  presented 
based  upon  gradually  varied  flow  theory. 
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2.  DEVELOPMENT  OF  THE  VERTICAL  STACK  TO  HORIZONTAL  DRAIN  NUMERICAL  MODEL 

2.i  ENTRY  ENERGY  BASED  ON  FULLY  DEVELOPED  VERTICAL  STACK  FLOW 

Figure  3 illustrates  the  annular  nature  of  vertical  stack  flow  under  fully 
developed  flow  conditions  [2]. 

For  a section  AL  of  the  stack  a force  balance  equation  for  the  falling  water 
may  be  expressed  under  steady  developed  flow  conditions,  where  stress  at  the 
air  water  interface  is  neglected,  as 

tiD  t0(AL)  = pir  Dt(AL)g  (1) 

where  D = pipe  diameter 

t = annular  thickness 

x0  = wall  shear  stress  resisting  motion 
It  follows  that: 


= P t g 


(2) 


Following  the  work  of  Colebrook  White  it  may  be  shown  that  the  shear  stress  may 
be  expressed,  for  full  bore  pipe  flow  by  the  equation 


/p~y2 

V S T0 


2 log10 


[ 


k 

3.7  D 


+ 


2.51 

Re 


(3) 


where  V is  flow  mean  velocity 

k is  surface  roughness  coefficient 
and  Re  is  Reynolds  Number,  based  on  pipe  diameter. 


Ackers  (3)  has  shown  that  equation  (3)  may  be  applied  to  partially  filled  pipe 
flow  if  the  characteristic  length  D is  replaced  by  the  hydraulic  radius: 


m 


Area 

Wetted  Perimeter 


which  for  a full  bore  pipe  flow  is 
m = D/4 

Replacing  Re  by  VD/ v and  substituting  for  D in  terms  of  m yields 


(4) 


P V^ 


8 x. 


= - 2 log10 


k + 2.51v  /p  V 
14.8  m 4 Vm  J 8 t 


o 


(5) 


For  the  annular  flow  condition  in  the  vertical  stack,  it  has  been  shown  that 
= P t g 


2 


and,  it  also  follows  that 


(6) 


if  t <<  D,  as  is  the  case  here  [2,4] 

Substituting  from  (2)  and  (6)  into  (5)  yields  an  expression  for  annular  thick- 
ness, t,  where  the  volumetric  flow  Q is  taken  as 


This  expression  may  be  solved  to  yield  t for  a given  Q and  hence  the  annular 
velocity  V. 

The  specific  kinetic  energy  in  the  annular  downflow  then  becomes 


Naturally,  the  expressions  above  only  apply  strictly  to  fully  developed  flow. 
Previous  work  [4]  has  shown  that  the  necessary  vertical  stack  length  to  reach 
this  condition  is  less  than  3m  in  many  design  cases,  i.e.,  less  than  a story 
height.  In  those  cases  the  derivation  above  may  therefore  be  utilized.  The 
influence  of  vertical  stack  pipe  wall  roughness  for  given  flow  rate  Q,  upon 
x 0 follows  from  (8)  and  (2). 

The  flow  rate  in  the  stack  may  vary  with  time  due  to  the  imposed  load  from  the 
branch  drain,  as  sketched,  figure  4.  However,  it  may  be  assumed  that  the  flow 
rate  profile  is  made  up  of  an  infinite  number  of  fully  developed  flow  sections, 
thus  allowing  the  annular  thickness  and  energy  vs.  time  profiles  to  be  devel- 
oped as  shown.  There  profiles  may  then  be  taken  as  the  entry  energy  profiles 
at  the  drain  upstream  boundary.  With  sufficient  height  of  fall  in  the  vertical 
stack  the  equilibrium  terminal  velocity,  Vt , is  attained. 

2.2  INITIAL  FLOW  PROFILES 

The  development  of  flow  depth  along  the  length  of  the  horizontal  drain  pipe 
under  the  initial  steady  flow  conditions  may  be  generated  by  reference  to  grad- 
ually varied  flow  theory.  The  basic  assumption  is  that  the  local  head  at  any 
section  is  given  by  the  equivalent  steady  flow  loss  equation,  either  Manning's 
or  Darcy's  equation,  utilizing  the  Colebrook  White  friction  factor  expression. 

Depth  profile  predictions  by  numerical  integration  are  based  on  this  assumption, 
expressed  in  figure  5 by 


Q = it  Dt  V 


(7) 


where  V is  the  mean  annular  velocity 


(8) 


2 g 


3 


d_ 

dL 


^ + (ZQ  - SoL>  + h 


= - S 


(9) 


where  S is  the  slope  of  the  energy  grade  line  determined  from  Mannings  equation 
or  from  the  general  Chezy  expression: 


v/T  = n Q 


S = _v 
C^m 


A m^/3 

2 


(Manning,  metric) 
(Chezy) 


(10) 

(11) 


Hence  - 1 dv  + ^ - dh  = S 

g dL  ° dL 

Then,  from,  Q = VA 

dV  A + V dA  = o 

dL  DL 

Also,  as  dA  = Tdh,  where  T is  the  water  surface  width  (which  does  not  change 
with  L,  or  changes  slowly)  it  follows  that 

dV  - _ V dA  - _ VT  dh  = _ QT  dh 

dL  A dL  A dL  A2  dL 

and  substituting  yields 


/ Q2T\  dh  + Sq 
VgAVdL 


- dh  = s 
dL 


4L  ■ i - s)  ,<-v] 


dh 


or  the  distance  along  the  flow  profile,  between  two  depths  hQ  and  hj , becomes 
hl  r~ 


- I 


1 ~ Q2  T/gA3 
(SQ-S) 


] 


dh 


(12) 


Figure  6 illustrates  this  integration. 

In  order  to  numerically  evaluate  (12)  it  is  necessary  to  define  a boundary 
depth  at  the  pipe  entry. 

The  available  flow  energy  at  the  horizontal  drain  entry,  assuming  no  losses,  is 
based  on  the  vertical  stack  annular  specific  flow  energy, 


E = h 


2 

la 

2g 


+ Xa 


(13) 
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where  hQ  and  VQ  are  the  flow  depth  and  velocity  at  entry  to  the  horizontal 
drain. 

1  * * * *  Vt2 

Substituting  E = — — — , yields 

2 g 


1 

2 


ho  + 


ho  + 


1_ 

2g 


= F(h0) 


(13a) 


and  solving  for  h0,  with  VQ  = Qq/Aq  and  A0  = f(hQ)  yields  the  upstream  bound- 
ary depth. 


The  downstream  boundary  may  be  taken  as  the  flow  normal  depth,  hn  at  the 
distance  downstream  to  be  determined  from  the  integration  of  (12).  As  the 
flow  is  supercritical,  the  depth  may  rise  gradually  to  this  value  with  no 
discontinuity,  such  as  a hydraulic  jump. 

The  choice  of  dh  values  in  the  numerical  integration  is  based  on  the  difference 
between  the  control  depth,  hQ  and  the  downstream  "target"  depth,  hn.  A reason- 
able integration  step  size  is  given  by 


dh  = (hn  - ho)/30 


(14) 


Once  the  depth  reaches  hn  the  remaining  downstream  profile  is  assumed  to  be  at 
that  constant  depth. 


For  short  pipe  lengths,  or  high  energy  steady  initial  flows,  it  may  be  that 
the  normal  depth  is  not  achieved  within  the  available  pipe  length.  For  this  case 
the  integration  is  terminated  when  SAL  exceeds  the  pipe  length,  i.e.,  ZAL>L. 

A valid  flow  profile  is  however  achieved  in  this  case. 


2.3  INITIAL  CONSIDERATION  OF  LOSS  COEFFICIENT  AT  ENTRY  TO  HORIZONTAL  DRAIN 


Figure  7 illustrates  the  conditions  for  an  initial  derivation  of  the  energy 

loss  at  the  entry  elbow  bend  to  the  horizontal  drain.  It  is  assumed  that 
within  the  annular  flow  the  specific  energy,  e,  along  "streamlines"  of  the 

downflow  is  represented  by  the  kinetic  energy  at  the  terminal  velocity  condi- 

tion. Degradation  of  the  energy  is  assumed  to  result  from  mixing  at  the  bend 
due  to  the  varying  orientation  of  flow  directions  from  the  vertical  to 
horizontal.  The  equation  may  be  assumed  as: 

V 2 

e = ~ (sin2  I o)  (14a) 

where  Vt  is  the  mean  annular  (terminal)  velocity.  The  equation  (14)  indicates 
the  vertically  falling  flow  at  0=0  loses  its  entire  energy  content  and  only 
adds  mass  at  the  drain  inlet  and  at  0=tt  the  streamline  total  energy  content  is 
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preserved.  Then  a relation  for  the  specific  energy  content  of  the  annular 
mass  flow*  may  be  obtained  as: 


/ e(|-  - t)  t d 0 = Efi  JL  £d2  - (D-2t) (15) 


Where  Eg  is  the  annular  energy.  Substituting  for  e yields 


£d2  - (D-2t)2J 

= 2 

1 (D-t)  t — f sin2  JLVQ 

0 2 2g  \ 2/ 

= (D  - t)  1 V2 
2 g C 

7 

0 

^sin2  ©JdO 

D t V’-2 

= (I  " g ' 

(tt/2) 

v 2 

„ (D/2-t)  Vt 

(D-t)  2g 


For  t « D,  Eg  + 

2 2g 


(16) 


The  model  does  not  apply  for  t -*■  D/2  (full  bore  pipe  flow)  since  the  turning 
losses  at  an  elbow  for  full  bore  flow  do  not  have  comparable  mixing.  Flow 
separation  from  full  bore  downflow  into  a partially  filled  horizontal  drain 
requires  different  physical  modeling. 


* The  term  (D/2-t)  includes  other  losses  for  annual  flow  not  in  contact  with 
the  pipe  wall.  The  use  of  (D-t) 1/2  (in  average  cross  section  area  formul- 
ation for  equation  (15))  leads  to  a constant  energy  loss  of  50  percent. 
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Values  of  the  factor  modifying  the  kinetic  energy  term,  Vtz/2g,  are  indicated 
below: 


Pipe  Diameter 
(mm) 


Annular  Film 
Thickness 


Ratio 

Factor 


(mm) 


(D/2-t)/(D-t) 


100 


25. 

12.5 

2.5 


33 

43 

49 


Measured  values  from  experiments  for  the  annular  thicknesses  are  difficult  to 
obtain  [2,  4];  however,  indicated  ranges  are  much  less  than  20  percent  of  the 
diameter  and  are  a function  of  flow  rate  in  the  pipe,  the  pipe  diameter,  and 
wall  roughness.  For  D>t,  (D/2-t)/(D-t)  -»■  0.5  hence  an  energy  "loss"  term  of 
50  percent  may  be  a reasonable  assumption  for  preliminary  studies.  Conditions 
other  than  for  terminal  velocity  should  retain  the  relationship  (16)  with 
replaced  by  the  falling  velocity  as  a function  of  length  along  the  stack. 

A study  of  entry  loss  effects  could  naturally  include  K values  greater  than 
unity.  That  condition  would  imply  that  the  entry  flow  energy  was  higher  than 
that  appropriate  to  gravity  driven  stack  flow  as  might  be  the  case  if  the 
drain  were  supplied  from  a pressurized  source.  Surge  flows  generated  in  tests 
reported  in  [5]  appear  to  correspond  to  pumped  head  conditions;  there,  other 
elbow  losses  may  also  be  significant  and  required  for  inclusion  of  the  analysis 
of  the  experimental  results. 

3.  FLOW  SIMULATION 

3.1  INITIAL  FLOW  DEPTH  PROFILES  ALONG  THE  HORIZONTAL  DRAIN 

Figures  8 and  9 illustrates  the  depth  profiles  along  a drain  set  at  slopes  of 
0.02  and  0.01  and  the  other  identified  parameters  for  a range  of  initial  flow 
rates.  The  profiles  attain  the  normal  depth  at  some  distance  along  the  pipe, 
thereafter  the  flow  is  assumed  to  continue  at  that  depth  until  discharge  at 
the  pipe  exit.  The  pipe  length  necessary  to  attain  the  steady  normal  flow 
depth  increases  for  larger  values  of  the  steady  flows  rates,  or  as  the  pipe 
slope  is  increased. 

Figures  10  and  11  illustrate  the  changes  in  profile  for  the  8 i/s  case  in  pipes 
pitched  at  1/100  and  1/50  for  both  a decreasing  and  increasing  inlet  energy 
factors.  The  reduction  in  entry  energy  indicated  by  the  introduction  of  a 
loss  coefficient  will  steepen  the  depth  profile  resulting  in  normal  depths 
being  achieved  earlier  in  the  pipe.  The  determination  of  the  correct  flow 
depth  profile  for  the  start  of  experiments  or  predictive  simulation  studies  Is 
important  since  it  effectively  provides  the  flow  depth  condition  immediately 
ahead  of  the  time-dependent  surge  wave. 

It  will  be  seen  that  if  the  inlet  energy  is  artificially  increased  (e.g.  pumped 
inlet  condition  [5]  then  the  normal  flow  depth  development  is  delayed.  Th< 
corollary  to  this  is  that  a fixed  length  pipe  may  be  seen  to  carry  rather  ;i  >r. 
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flow  than  suggested  by  the  steady  flow  calculations  of  normal  depth.  For 
example,  with  K > 1.0  with  pipe  lengths  less  than  16ra,  the  depths  are  less 
than  those  which  occur  with  gravity  drain  conditions,  from  which  it  would  be 
implied  that  higher  flow  rates  can  exist  in  the  pipe  without  violating  permis- 
sible half-filled  depths.  Care  should  be  taken  in  any  experimental  work  in 
this  area  to  ensure  that  the  entry  energy  is  appropriate.  Failure  to  do  this 
could  suggest  enhanced  carrying  capacities  that  will  not  be  found  in  practice. 
The  presence  of  pipe  fittings,  such  as  a junction  will  result  in  substantial 
energy  dissipation  and  the  flow  will  therefore  rapidly  increase  in  depth, 
could  fill  the  pipe  cross  section  and  be  improperly  interpreted  as  the  hydraulic 
jump. 

Table  1 presents  normal  depth  values  for  100  mm  diameter  pipes  at  gradients  of 
1/100  and  1/50. 

Table  1.  Normal  Depths  in  a 100  mm  Diameter  Smooth  Drain, 
k = 0.0,  Calculated  via  the  Colebrook  White 
Friction  Loss  Expression 


Pipe  Gradient  Flow  Rate 

SL/s 


Normal  Depth 
mm 


1/100 


1.0 

23.6 

2.0 

33.6 

3.0 

41.8 

4.0 

49.2 

5.0 

56.2 

6.0 

63.2 

7.0 

70.8 

8.0 

79.4 

1.0 

19.5 

2.0 

27.5 

3.0 

33.9 

4.0 

39.5 

5.0 

44.8 

6.0 

49.7 

7.0 

54.5 

8.0 

59.3 

9.0 

64.1 

10.0 

69.2 

11.0 

74.7 

12.0 

81.1 

13.0 

91.7 
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3.2  SIMULATION  OF  VERTICAL  STACK  ANNULAR  THICKNESS  AND  TERMINAL  VELOCITY 

Equation  8 may  be  utilized  to  calculate  the  annular  flow  thickness  in  the 
vertical  stack  for  a range  of  steady  flow  rates.  The  terminal  velocity 
corresponding  to  this  thickness  and  flow  rate  may  then  be  calculated  from 
equation  7. 

Figures  12  and  13  illustrates  the  variation  in  terminal  velocity  and  annular 
thickness  with  flow  rate  and  also  include  previous  experimental  work  reported 
elsewhere  by  Wise  and  Wyly  [2,4].  The  experimental  values  of  terminal  velocity 
probably  exceed  those  predicted  due  to  the  difficutly  in  measuring  the  thick- 
ness of  the  annular  layer  (the  air  core  interface  with  the  water  edge  is  at  the 
same  velocity  due  to  the  shear  interaction)  and  as  some  of  the  flow  will  be 
’’lost"  as  spray  falling  down  the  center  of  the  pipe.  Thus  the  value  of  terminal 
velocity  calculated  from  equation  7 will  be  excessive  as  t is  underestimated. 
Wise's  results  also  indicate  this  tendency. 

With  no  energy  loss  in  the  4'ntry  bend  the  entry  pipe  flow  at  velocity  is  close 
to  the  stack  terminal  velocity.  The  results  presented  in  figures  13  and  14 
conform  to  the  energy  boundary  equation  (13a). 

3.3  SURGE  ATTENUATION  FOLLOWING  STACK  DISCHARGE  TO  A HORIZONTAL  DRAIN 

The  necessary  equations  to  simulate  wave  attenuation  along  the  length  of  the 
pipe  have  been  developed  [6] . The  transition  from  the  stack  to  the  horizontal 
drain  may  be  characterized  by  the  entry  energy  boundary  condition  described 
previously,  based  on  the  annular  terminal  flow  in  the  stack.  Initial  condi- 
tions along  the  horizontal  drain  may  be  set  up  utilizing  the  gradually  varied 
flow  theory  outlined. 

The  wave  formation  in  the  horizontal  drain  may  be  viewed  as  being  made  up  of 
two  separate  but  superimposed  mechanisms.  As  the  surge  progresses  down  the 
drain  its  depth  increases  as  energy  is  dissipated.  At  the  same  time  the  wave 
attenuates  due  to  the  diffusion  so  that  the  superposition  of  these  two  effects 
results  in  a peak  flow  depth  downstream  of  the  entry  along  the  pipe  length. 
Figures  15  and  16  indicate  the  peak  flow  depths  predicted  for  a range  of  inflow 
profiles.  Figures  17  and  18  present  the  calculated  peak  flow  rates,  calcu- 
lated at  each  pipe  section  as  the  product  of  (VxA).  It  will  be  seen  that 
the  peak  flow  rates  exceed  the  maximum  inflow  rates.  This  apparent  anomaly  is 
explained  as  the  rising  inflow  surge  generates  waves  moving  with  increasing 
wave  speed  greater  than  the  mean  velocity.  The  waves  "overtake"  the  preceding 
flow  and  hence  generate  both  depths  and  flow  rates  in  excess  of  those  predicted 
by  a steady  flow  analysis. 

It  will  also  be  seen  that  any  surge  alleviation  provided  by  the  reduction  in 
surge  flow  is  only  manifested  for  the  triangular  inflow  profiles.  Profiles 
displaying  a plateau  at  peak  inflow  rate  tend  to  produce  nearly  similar  maxi- 
mum depths.  This  would  suggest  that  the  increased  flow  capacity  often  referred 
to  as  a result  of  attenuating  surge  flows  is  applicable  only  to  relatively 
"sharp",  triangular  inflow  rate  vs.  time  profiles. 
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Figures  19  and  20  indicate  the  location  of  the  peak  flow  and  peak  depth  along 
the  horizontal  drain  for  a range  of  inflow  profiles  (the  calculation  nodes  are 
approximately  0.5m  apart).  It  will  be  noted  that  the  peak  flow  rate  invari- 
ably preceeds  the  location  of  peak  depth.  The  same  effect  is  also  noticeable 
in  figures  21  and  22  that  indicate  the  peak  depth  and  flow  rate  "envelopes" 
along  the  horizontal  drain  at  gradients  of  1/50  and  1/100.  Figures  23  and  24 
illustrate  the  flow  depths  at  three  stations  along  the  drain  over  the  duration 
of  the  simulation. 

3.4  EFFECT  OF  BASEFLOW  ON  SURGE  ATTENUATION 


Figures  25  and  26  present  peak  flow  depths  and  rates,  together  with  their 
location  for  two  computational  test  cases,  for  baseflows  0.2  £/s  and  2 Si/s 
with  triangular  flow  profiles  rising  to  12  1/ s. 

It  would  appear  from  these  results  that  the  baseflow  only  becomes  an  important 
parameter  when  the  inflow  surge  time  is  sufficiently  long  so  that  the  benefi- 
cial attenuation  of  the  peak  does  not  occur.  Reference  to  the  earlier  figures 
15  and  16  indicate  that  only  the  2 second  base  triangle  surge  suffered  any 
appreciable  attenuation  below  the  steady  flow  normal  depth  associated  with 
its  peak. 

It  is  apparent  that  further  data  is  required  to  be  more  specific  as  regards 
the  effects  of  baseflow. 

3.5  COMPARISON  OF  PREDICTED  AND  OBSERVED  DEPTH  PROFILES  FOR  A HORIZONTAL  DRAIN 
SUPPLIED  BY  A VERTICAL  STACK 

Figures  27  and  28  present  comparisons  between  experimentally  measured  depth  vs. 
time  profiles  in  a 75  mm  pipe  set  at  a gradient  of  0.06.  Surges  enter  the  pipe 
via  a vertical  stack  supplied  by  a constant  baseflow  plus  discharges  from  water 
closets  at  higher  levels.  The  results  presented  were  recorded  in  a dynamics 
of  plumbing  system  test  program  [7,8]. 

Unfortunately  no  data  is  available  on  the  soil  stack  or  drain  inflow  entry 
flow  rate  profile  which  is  governed  by  the  mixing  at  the  horizontal  drain 
branch  elbows  into  the  soil  stack  downflow.  The  flow  rate  profiles  presented 
in  figures  27  and  28,  for  the  two  baseflow  cases,  were  best  estimates  of  the 
likely  water  closet  discharge  profiles  determined  from  the  measured  depth  vs. 
time  records  in  the  horizontal  branch  drains  at  the  second  and  third  (story) 
level  leading  to  the  vertical  stack.  A single  water  closet  on  the  third  level 
and  two  back-to-back  water  closets  on  the  second  level  were  simultaneously  dis- 
charged into  constant  downflows  in  the  vertical  stack. 

The  estimated  discharge  profiles  were  summed  (as  shown  in  the  figures)  and 
used  to  establish  the  inlet  conditions  to  the  drain  model  based  upon  the  ter- 
minal velocity  kinetic  energy  equation.  Comparison  between  observed  and 
predicted  depth  variation  at  a downstream  depth  tapping  are  shown.  The  time 
scale  of  the  observed  results  have  been  adjusted  so  that  the  observed  and 
predicted  profiles  move  from  the  steady  base  flow  condition  at  the  same  time. 
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This  approximate  adjustment  may  account  for  some  of  the  time  scale  deviations 
shown.  The  predicted  normal  depths  with  the  steady  baseflow  (after  the  surge) 
are  within  2 percent  of  measured  values. 

The  results  presented  in  figures  27  and  28  support  the  proposed  model  for 
stack  to  horizontal  drain  entry  model.  The  baseflow  does  not  affect  the 
predicted  peaks  appreciably;  however  that  was  expected  from  the  results  shown 
in  figures  25  and  26.  Any  improvement  on  the  agreement  requires  an  improved 
knowledge  of  the  inflow  flow  rate  profile,  or  at  least  a better  estimate  of 
the  water  closet  flow  rate  discharge  curves  as  modified  by  the  junction  fit- 
tings at  the  upper  levels  from  the  horizontal  branch  drains.  Further  research 
is  necessary  to  fully  account  for  mixing  at  entry  fittings  into  the  soil  stack 
downflow.  There  have  been  no  measurements  of  the  distribution  of  annular 
flow  properties  in  the  vertical  soil  stack  which  are  required  to  provide 
understanding  of  such  gravity  driven  flow  regimes.  Investigations  of  annular 
flow  phenomena  indicate  that  surface  waves  are  observed  on  vertical  falling 
films.  That  suggests  that  superposition  of  surges  on  steady  soil  stack  flows 
may  take  place,  wave  velocities  greater  than  the  mean  falling  flow  velocities 
can  occur,  and  flow  redistribution  may  result,  in  addition  to  mass  mixing  at 
junctions  or  fittings. 

For  the  same  surge  profile  conditions  of  figure  27,  the  predicted  effects  of 
the  entry  loss  factor  for  K=0.5  and  1.0  are  shown  in  figure  29  at  three  sta- 
tions along  the  length  of  the  drain.  The  higher  inlet  energy  condition 
(K=1.0)  result  in  smaller  peak  depths  compared  with  the  lower  inlet  energy 
condition  (K=0.5).  A steepened  surge  wave  develops  with  higher  energy  content, 
due  to  the  higher  wave  velocity  (/gA/T)  overtaking  the  upstream  flow  which 
reinforces  the  surge  surge  wave  front.  Also  it  is  noted  that  the  surge  depth 
attenuation  for  the  higher  energy  case  has  not  yet  occurred  up  to  the  10  m 
station.  At  that  station  the  1/ D ratio  (length  of  pipe  to  pipe  diameter)  is 
about  130;  the  prediction  of  large  values  of  2,/D  were  reported  in  earlier 
experiments,  reported  in  [5]  with  pumped  inlet  conditions.  The  sharp  falloff 
from  peak  depth  for  the  high  energy  computations  at  the  10.0  m station  indi- 
cates that  higher  wave  velocities  have  overtaken  the  surge.  Those  comparisons 
indicate  that  a need  exists  to  improve  the  simplified  drain  inlet  energy  profile 
derivation.  Observations  made  in  experiments  [6,  7]  suggest  that  a transition 
entry  length  in  the  horizontal  drain,  downstream  of  the  elbow,  may  be  required 
to  account  for  a flow  mixing  transition  in  the  modeling. 
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In  figure  30,  the  results  for  the  drain  depths  for  three  inlet  profiles  are 
shown.  The  inlet  profiles  considered  are  indicated  below. 


Flow 


The  assumed  conditions  represent  waste  discharge  loading  patterns  due  to  a 
short  time  interval  surge  followed,  after  an  interval,  by  a rise  to  a steady 
discharge  flow.  In  all  cases,  the  initial  surge  rapidly  attenuates  (as  has 
been  noted  for  figures  15  and  16).  With  increasing  times  the  transient  depths 
occur  greater  than  the  steady-state  conditions  before  trending  towards  the 
normal  depths  (compare  with  values  of  table  1).  For  case  A,  the  peak  depths 
are  90  percent  of  full  bore  flow,  whereas  the  normal  depth  was  shown  to  be 
about  80  percent  of  full  bore.  Currently  accepted  plumbing  code  criteria  for 
allowable  partially  filled  drain  pipe  depths  permit  only  half-filled  pipes 
(although  the  report  [9]  on  which  the  codes  are  based,  assumed  75  percent  for 
surge  flows).  Neither  the  transient  measurement  methods  or  calculation  methods 
were  available  to  the  investigators  or  code  writing  panels  at  that  time, 
consequently  the  effects  discussed  herein  were  only  available  by  observations  and 
greatly  simplified  measurement  methods.  Here,  the  larger  wave  speeds  at  the 
greater  rate  of  increase  to  the  higher  steady  flow  rates  cause  a steepening  of 
the  wave  and  occurrence  at  somewhat  earlier  times  for  the  same  station.  The 
peak  depths  occur  at  pipe  length  to  diameter  (A/D)  ratios  of  about  100  for  case 
A,  and  75  for  cases  B and  C.  The  predicted  results  are  consistent  with  the 
data  from  the  pumped  experiments  [5]  that  showed  the  limit  of  hydraulic  drain 
capacity,  based  upon  a prescribed  depth  for  a fixed  length,  is  limited  by  lowest 
steady  flow  rate  surge  component  and  similar  ranges  of  (A/D)  values. 
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4.  CONCLUSIONS 


This  report  has  introduced  the  concept  of  an  energy  entry  condition  to  simulate 
vertical  stack  flow  entering  a horizontal  drain.  In  general  the  model  has  been 
satisfactory  in  that  no  numerical  instability  problems  have  been  encountered. 
The  predictions  of  flow  depth  and  flow  rate  through  the  stack  and  drain  match 
the  physical  mechanisms  known  to  occur  for  these  flow  situations.  The  predic- 
tive model  provides  for  the  first  time,  a comprehensive  theoretical  basis  for 
explanation  of  the  earlier  experiments  at  NBS  to  establish  the  hydraulic  capa- 
city of  drains. 

It  may  be  concluded  that  the  simplified  annular  stack  flow  and  entry  energy 
model  provides  the  first  theoretical  basis  as  a suitable  method  of  linking 
networks  of  a number  of  horizontal  drains  with  vertical  stack  discharging  to 
a lower  level  pipe.  Network  horizontal  pipe  junctions  solutions  have 
progressed,  in  a parallel  study  [10],  so  that  complex  interconnected  drainage 
system  analysis  is  feasible  now  through  combination  of  the  predictive  models. 
The  models  will  be  developed  further  and  validated  during  the  remainder  of  the 
projected  program. 
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Figure  1.  Summary  of  method  of  characteristics  solution  of  the  unsteady 
partially  filled  pipe  flow  equations 
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Figure  2.  Application  of  the  method  of  characteristics  to  unsteady  super- 
critical partially  filled  pipeflow 
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Figure  3.  Fully  developed  annular  flow  in  vertical  stacks 
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Figure  4. 
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Figure  5.  Calculation  of  the  energy  grade  line  for  gradually  varied  flow 
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Figure  6.  Development  of  flow  depth  profiles  in  gradually  vared  pip  flow 
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Figure  7c  Upstream  boundary  formed  by  an  entry  energy  equation 
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Figure  8.  Depth  profiles  for  an  initial  steady  flow  along  the  drain  at  a 
1/50  gradient,  no  energy  loss  at  entry 
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Figure  9.  Depth  profiles  for  an  initial  steady  flow  along  the  drain  at  a 
1/100  gradient,  no  energy  loss  at  entry 
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Figure  10.  Influence  of  the  entry  energy  loss  coefficient  on  initial  steady 
flow  depth  profiles  (1/100) 


Figure  11.  Influence  of  the  entry  energy  loss  coefficient  on  initial  steady 
flow  depth  profiles  (1/50) 
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Annular  stack  flow  thickness  as  a function  of  stack  steady  flow  rate 


Figure  13.  Stack  annular  flow  terminal  velocities  compared  to  existing 
experimental  data 
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Figure  14.  Flow  depth  at  pipe  entry  compared  to  annular  stack  flow  tM 
for  a range  of  steady  flow  rates 
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Figure  15.  Peak  depths  predicted  along  the  horizontal  drain  for  a range  of 
inflow  profiles  (1/50) 
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Figure  16.  Peak  depths  predicted  along  the  horizontal  drain  for  a range  of 
inflow  profiles  (1/100) 
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Figure  17.  Peak  Flow  rates  predicted  locally  along  the  horizontal  drain  for  a 
range  of  inflow  profiles  (1/50) 
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Figure  18.  Peak  flow  rates  predicted  locally  along  the  horizontal  drain  for  .1 
range  of  inflow  profiles  (1/100) 
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Figure  19.  Location  of  predicted  peak  flow  depth  and  local  flow  rate  for  a 
range  of  inflow  profiles  ( 1 / 50 ) 
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Figure  20.  Location  of  predicted  peak  flow  depth  and  local  flow  rate  for  a 
range  of  inflow  profiles  (1/100) 
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Figure  21.  Peak  flow  depth  and  local  flow  rate  envelopes  along  the  length  of  the  horizontal  drain  ( 1 / SO ) 
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Figure  ??.  Peak  flow  depth  and  local  flow  rate  envelopes  along  the  length  of  the  horizontal  drain  (1/1 00 ) 
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Figure  23.  Depth  vs.  time  profiles  at  a number  of  stations  along  the  horizontal  drain  (1/50) 
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number  of  stations  along  the  horizontal  drain  (1/100) 
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Figure  25.  Summary  of  baseflow  effect  on  peak  depth  and  local  flow  rate  predictions  (1/50) 
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Fleur*  2A.  Summary  of  haseflow  effect  on  peak  depth  end  local  flow  rate  predictions  (I /inn) 


Figure  27.  Comparison  of  observed  and  predicted  surge  attenuation  downstream 
of  a vertical  stack  to  horizontal  drain  entry 
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Figure  28.  Comparison  of  observed  and  predicted  surge  attenuation  downstream 
of  a vertical  stack  to  horizontal  drain  entry 
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Figure  29.  Comparison  of  predicted  surge  depths  downstream  of  a vertical  stack 
to  a horizontal  drain  entry  for  two  entry  energy  loss  factors 
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Figure  30.  Comparison  of  predicted  surge  profile  depths  in  a horizontal  Ir.iln 
downstrem  of  a vertical  stack  for  composite  loading  profiles 
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